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5^ . Abstract 

We analyze the production and subsequent decay of the neutral Higgs bosons h = h", H", A" 
I of the MSSM into electrically neutral quark pairs of different flavors {qq' = tc, bs, depending 

^ ■ on h) at the LHC, i.e. <j{pp — > /i ^ qq'), and compare with the direct FCNC production 

CO ' mechanisms cr(pp qq')- The cross-sections are computed in the unconstrained MSSM with 

minimal flavor- mixing sources and taking into account the stringent bounds from b — > S7. We 
extend the results previously found for these FCNC processes, which are singularly uncommon 
in the SM. Specifically, we report here on the SUSY-EW part of cr(pp h ^ qq') and the 
|0 I SUSY-QCD and SUSY-EW contributions to a{pp — > bs). In this way, the complete map of 

■ MSSM predictions for the gg'-pairs produced at the LHC becomes available. The upshot is 
I that the most favorable channels are: I) the Higgs boson FCNC decays into bs, and 2) the 

■ direct production of tc pairs, both of them at the ^ Ipb level and mediated by SUSY-QCD 
effects. If, however, the SUSY-QCD part is suppressed, we find a small SUSY-EW yield for 
(t(pP h te)max ^ lO^^pb but, at the same time, (T(pp — > /i — s- fes)max ^ O.I — 1 pb, which 



o 



X 



M ' implies a significant number 10* — 10^) of bs pairs per 100 fb ^ of integrated luminosity. 



1 Introduction 



The upcoming generation of TeV-class colliders, headed by the imminent startup of the Large 
Hadron Collider (LHC) at CERN, will offer us the opportunity to dig deeper than ever into the 
nature of fundamental interactions. Despite its successful career, the Standard Model (SM) of 
elementary particles still offers a number of intriguing puzzles to be resolved, such as the ultimate 
origin of the electroweak (EW) symmetry breaking and the mass generation mechanism. Among 
the wide set of theoretical proposals that have been conjectured so far, the possibility that the 
fundamental laws of Nature exhibit a symmetry between fermionic and bosonic degrees of freedom 
- that is, supersymmetry (SUSY) and most particularly the Minimal Supersymmetric Standard 
Model (MSSM) [1] - has been postulated as one of the firmest candidates to account for possible 
scenarios of physics beyond the SM, i.e. of New Physics (NP). 

The quest for experimental signatures of an underlying SUSY dynamics is undoubtedly one of 
the major endeavors for the LHC. Besides the direct production and tagging of SUSY particles 
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(which could be cumbersome in practice), it is of great interest to consider the impact of SUSY 
radiative corrections on conventional processes in the SM. In such cases the enhancement capabili- 
ties associated to the non-standard couplings may provide sources of genuine NP signatures. This 
possibility was thoroughly analyzed in the past for the W and Z boson physics (sec e.g. [2,3]). 

More recently, a great deal of attention has been devoted to processes triggered by Flavor- 
Changing Neutral-Current (FCNC) sources beyond the SM, in particular from SUSY interactions. 
From the seminal work of Glashow, Iliopoulos and Maiani (GIM) [4] it has been known that the 
FCNC interactions are absent at the tree level within the SM and, most significantly, they are 
largely suppressed at the 1-loop order. This, so-called, GIM mechanism is an in-built feature of 
the mathematical core of the SM, thanks to the unitarity of the CKM matrix. Low energy B-meson 
physics, for example, provides a number of well-measured FCNC processes, such as the B-meson 
radiative decay b — > 37. Its branching ratio reads B{h — > S7) = (2.1 — 4.5) x 10~^ at the So- 
level [5]. This important and well studied process can be used to constrain models of NP that 
predict non-standard flavor-changing interactions. In contrast, the FCNC effects involving the top 
quark as an external particle turn out to be dramatically suppressed by the GIM mechanism. The 
predicted branching ratios are at the level of B{t — eg) ~ 10"-*^^ and B{t cH) ~ lO"-'^^ [6], 
hence far below the limits of observability. 

In stark contrast with this meager SM panorama, the MSSM opens new vistas for a fruitful 
FCNC physics of the top quark. A most remarkable example is the case of the top quark decay 
into the lightest supersymmetric neutral Higgs boson, h°, whose branching ratio could be of order 
B{t ^ ch°) ~ 10~^ — 10~^ [7] - previously underestimated in [8]. A branching ratio of this order 
represents not only an enhancement of 10 orders of magnitude above the SM prediction, but it 
opens the realistic possibility for measurement. The origin of these possible effects lies in the richer 
diversity of sources of flavor mixing in the MSSM, in particular among fermion and sfermions of 
the same charge and different generation. They ultimately stem from the so-called misalignment 
of the squark mass matrices with respect to the quark ones and can be described by means of 
the parameters S^^ (being A,B = L,R the chirality indices and i,j = 1,2,3 the flavor ones). 
A generic entry of the soft SUSY-breaking squark mass matrix reads (M^)4-S = S^^ fhi fhj (for 
i ^ j). As a consequence, new types of FCNC couplings arise and they need not be subdued by 
GIM suppression. Obviously, such scenarios bring us some hope to effectively unearth hints of 
SUSY physics out of the study of both low energy and high energy FCNC processes [9]. 

Interestingly enough, not only SUSY can help here; other alternative proposals of extended 
physics beyond the SM, among them the general Two-Higgs-Doublet Models (2HDM) [10], top- 
color models and strong flavor-changing effects predict in some cases an enhanced, and sometimes 
distinctive, FCNC phenomenology [11-14]. Attentive studies on this field can thus be of great help 
in seeking for new signatures and eventually disentangling the sort of NP hiding right there. ^ 

In this letter, we wish to further elaborate on the rich phenomenology of flavor-changing pro- 
cesses in the MSSM, specifically on those that could lead to an unsuspected overproduction rate of 
electrically-neutral heavy-quark pairs qq' of different flavors at the LHC. Most conspicuously, we 
have the process of single top quark production induced by FCNC, i.e. any of the two-body final 
states pp — ^ tc = tc -I- t c. This process is extremely suppressed in the SM, while it can be highly 
enhanced in the MSSM [16-20]. A similar situation applies to producing FCNC single b-events, 
i.e. bs = bs + sb, although here the SM suppression is not so drastic as in the tc case. 

There are different supersymmetric mechanisms leading to enhanced tc and bs final states. As 
a first possibility, we have the SUSY flavor-changing charged (and neutral) currents contributing to 
the direct production processes pp tc and pp ^ bs. An alternative mechanism to produce these 
final states is through the production and subsequent FCNC decay of a neutral Higgs boson, i.e. 

^For a review of top quark and Higgs boson FCNC physics in the MSSM and in general 2HDM models, see 
e.g. [15] and references therein. 
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pp — > /i ^ tc and pp ^ h ^ bs. This mechanism was explored within the context of the general 
2HDM [12], in which h = h^, H^, AP [10], and it was later revisited within the unconstrained MSSM 
through the strong supersymmetric (SUSY-QCD) corrections to the Higgs boson FCNC branching 
ratios into the heavy-quark final states [21,22]. However, a first computation of the corresponding 
LHC production rates appeared only in [23] . 

Supersymmetric electroweak (SUSY-EW) contributions to neutral Higgs boson decays could 
be important too. These effects can be induced, in principle, by both charged and neutral currents. 
The latter are triggered by the neutralinos and some authors have dealt with these effects [24-26] . 
Here, however, following [7] (see also [27]), we compute the SUSY-EW contributions in the super- 
CKM basis, i.e. under the assumption of minimal flavor violation (5^^ =0, i 7^ j). We argue 
and verify that the neutralino effects are subleading. Therefore, we concentrate on the charged 
current sector of the electroweak MSSM Lagrangian [28]. In this setup, we perform a systematic 
analysis along the lines of [23] by addressing the SUSY-EW effects on a{pp h ^ qq') induced 
by charginos, squarks and charged Higgs bosons in the MSSM. Wc also compute the pay-off in 
the number of hs events from the direct production mechanism (T(pp — > hs) within the MSSM, 
including both SUSY-QCD and SUSY-EW effects. Notice that the interest of knowing the SUSY- 
EW effects in these studies is that they could be the only sizeable supersymmetric source of FCNC 
gg'-pairs at the LHC, if gluinos turn out to be very heavy and/or the inter-generational mixing 
parameters would be too small or simply zero^. 

2 Higgs boson decays into neutral heavy quark-pairs 

In this section, we discuss the general expectations on the production cross-section of electrically 
neutral pairs of heavy quarks of different flavors at the LHC, whose origin stems from the FCNC 
decays of a neutral MSSM Higgs boson [10]. The main contribution to the overall process may be 
sequentially split into: i) production of a real neutral Higgs boson in a proton-proton collision, ii) 
followed by its decay through loop diagrams contributing to the FCNC final state. Assuming that 
the Higgs bosons are produced on-shell, the total production rate can be factorized as follows: 



Here V{h qq') is the total FCNC two-body partial decay width of the corresponding MSSM Higgs 
boson h = H*^, into the (kinematically possible) neutral states 6s = bs + sb or tc = tc -|- tc; 
and T{h ^Y) stands for the - consistently computed - total decay width in each case. 

Before coming to grips with the full numerical analysis, it is always a good exercise to perform 
an analytical estimate based on physical considerations. In doing this, we expect to gain insight 
into the dominant FCNC sources of enhancement within the MSSM and roughly predict (within 
order of magnitude, hopefully) the relation between the cross-sections a{pp h ^ tc) and (T(pp — > 
h ^ hs). In order to address such estimate, we first need to specify the possible supersymmetric 
electroweak contributions to the FCNC Higgs boson decays into the heavy-quark final states under 
consideration. For the sake of definiteness, let us concentrate on H*^ (i.e. the heavy CP-even Higgs 
boson of the MSSM). A sample of Feynman diagrams from the SUSY-EW sector is displayed in 
Fig. 1. The overall set is, of course, UV finite (as we have checked). However, we note that 
the subsets of charged gauge and Higgs boson mediated diagrams, on the one hand, and the 

^When considering the explicit FCNC effects on our processes, we just focus on 823 = 823 ■ We assume that the 
FCNC mixing may hold in the chiral LL sector of the squark 6x6 mass matrices only, which is well motivated 
theoretically. The other chirality sectors give similar contributions [7]. See also [16,18]. 



(7(pp — > /t — ^ qq') = cr(pp — ^ hX)B{h — > qq') 



= a{vv^hX)^ 




{qq' = bs or tc) . 



(1) 
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Figure 1: Sample of SUSY-EW contributions to the decay of the Higgs boson into bs from the 
charged current. The electroweak effects that could originate from the neutralinos are not shown. 
A similar collection of diagrams describes the corresponding decay into tc. 



chargino mediated, on the other, are separately UV-finite. Moreover, each of these subsets involves 
two different topologies, vertex corrections (V) and wave-function (WF) renormalization, whose 
respective contributions arc also separately UV-finite after summation over generations. 

It turns out that the chargino-mediated pieces drive the bulk of the contribution, basically 
because they are not suppressed by the GIM mechanism. Let us first consider the WF correc- 
tions. Using no other tools than dimensional analysis, power counting, CKM matrix elements and 
dynamical features, we can estimate the (one-loop) decay width of H° bs as follows: 

and similarly for the H*^ tc channel: 

Let us briefly explain the origin of the different terms appearing in the expressions above. In both 

cases we include the standard numerical factor arising from the loop integration. Gp stands for 
the Fermi constant, and Vtb ^ ^ 0.02, ^ 0.02 are the CKM matrix elements; MsuSY 

defines a common scale for the soft SUSY-breaking masses (of squarks and gauginos), the MSSM 
trilinear couplings are labelled by A^, At; and ji indicates the higgsino mass parameter. The latter 
accounts for the helicity flip that appears along the chargino linc^. The angle a is the mixing angle 
between the neutral CP-even Higgs boson states. Finally, tan/3 = V2/V1 is the ratio of vacuum 
expectation values of the two MSSM Higgs doublets [10]. A key element in the expressions above 
is the factor associated to the chiral transition of top squarks, mt (^t — tan /3) , or bottom squarks, 
mfe (^5 — yutan/?), depending on the case. In the estimate, we also include the Higgs-quark-quark 
and (the leading part of) the chargino-quark-squark couplings [10], from which we get the remaining 
dependences in a, 13. Last but not least, we also keep track of the relevant mass scales arising from 
the 2-point loop functions and the phase space integration, which we settle in terms of MguSY 
and the mass of the decaying Higgs particle, niH- Note that equation (2), for example, patently 
reveals different sources of possible large enhancements, mainly associated with large values of 

^For an explanation and a computation of the helicity flips in the quark propagator, see e.g. Refs. [29], speciflcally 
Fig. 7 and eq. (76) of the first work in [29]. 



4 



the trilinear coupling A^, the higgsino mass parameter fi and/or the trigonometric factor 1/ cos^/? 
(which behaves as ~ tan'^/3 in the large tan/3 regime). 

Similar arguments can be applied to the vertex corrections, and we arrive at 



3 (VtbVC^y ( mtml \ 



2 ^. \ ^ (■^^ sin a — /xcos a)"' ^ 
sin^ /? cos /? 

jjL {A}, cos a — jjL sin a) 
sin /3 cos^ /3 



(1) 

2 

(5) 



Some differences between the formulae above corresponding to tc and 6s final states are worth 
noticing. For instance, the quark (or squark) mass insertions involve distinct mass factors: mf, (in 
the tc case) and mt (in the 5s one). The trigonometric couplings are also different and for this 
reason one of the channels may be much more suppressed than the other at different regimes. In 
particular, we sec that in both cases the decay rate increases with tan /?, but while the leading 
effect of the bs channel lies in the WF corrections the dominant one in the tc channel resides in the 
V contribution. As for the neutralino yield (whose diagrammatic effects are not shown in Fig. 1), 
similar analytical estimates can be produced. We limit ourselves to single out some differences. 
To start with, we emphasize that these contributions arc proportional to the inter-generational 
mixing parameters and are, therefore, vanishing for 5f-^ = (i 7^ j). The corresponding effects in 
the WF sector are such that e.g. the factors Ai, — jj, tan /3 and At — 11/ tan (5 become interchanged in 
the counterparts of Eqs. (2)- (3) . Similarly with the factors At sin a — cos a and A}, cos a — jJLSma 
in the corresponding vertex contributions, i.e., the analogous of Eqs. (4)- (5). 

In sections 4-5, we will come back to the previous analytical estimates and shall compare them 
with the exact numerical results. This will be useful to track the dynamical origin of the leading 
MSSM contributions and, in particular, to argue that the neutralino effects are negligible. 

Similar considerations can be done for the other Higgs bosons. However, we recall that h'^ (the 
lightest CP-even Higgs boson) cannot decay into tc because m^o < mt in the MSSM. Moreover, 
for this Higgs boson, the implications of the so-called "small Oe//" scenario (triggered by large 
radiative corrections in the parameters of the MSSM Higgs sector) must be taken into account [30]. 



3 Framework for the numerical analysis 

In order to compute the SUSY-EW one-loop diagrams contributing to the FCNC cross-sections 
'^(PP — ^ — qq') for the three MSSM neutral Higgs bosons {h = h'^,H'',A°), we shall closely 
follow the notation and methods of Refs. [7,11,12,22,23]. Wc address the reader to these ref- 
erences for the technical details. In particular, a thorough exposition of the relevant interaction 
Lagrangians and similar set of Feynman diagrams for the FCNC interactions, is provided in [7]. 
See also Ref. [10] for basic definitions in the MSSM framework and [29] for detailed computational 
techniques and further illustration of the super symmetric enhancement effects in other relevant 
Higgs boson processes. Along our computation we have made use of HIGLU, PPHTT [31], LoopTools, 
FeynArts and FormCalc [32]. 

For the numerical evaluation, we shall adhere (wherever possible) to the general procedure put 
forward in Ref. [23], where significant parts of the numerical contributions to the cross-sections 
have already been reported. Next we specify the framework under which the computation of the 
observables (1) has been carried out in the present work: 

• We compute the one-loop SUSY-EW contributions to the FCNC partial decay widths T[h — ^ 
qq') and compare with the corresponding SUSY-QCD effects [23]; 
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• The SUSY-EW part is defined to be the set of contributions from charginos, neutrahnos and 
Higgs bosons. In practice, however, we will provide the detailed results from the charged 
current effects, and argue (and numerically verify) that the neutralino contributions are 
comparatively negligible. 

Our major goal is to assess what are the theoretical expectations on the observables (1) within 
the MSSM and, most particularly for the present work, to ascertain whether in the absence of 
strong supersymmetric sources of FCNC the electroweak SUSY sector can still provide significant 
rates. To this aim we have performed a systematic scan all over the MSSM parameter space and 
have determined the maximum values for the production rates (1) under study. Furthermore, in 
order to keep the CPU-time under a feasible range, such a scan has been carried out by means 
of a Monte Carlo sampling method [33] based on the well-known Vegas integration program [34]. 
This numerical procedure was amply tested in the similar computation presented in Ref. [23]. 

On this basis, we have performed a maximization of the FCNC cross-section within the following 
restrictions: 



tan/3 


50 


At 


\At\< 3MsusY 




\Ab\ < 3MsusY 




(0 • • • 1000) GeV 




"^Jl - ^dn - -^~9- -^SUSY 




-^SUSY 


MsvSY 


(150 • • • 1000) GeV 


nijifi 


(100 • • • 1000) GeV 


Ma- 

HI 


2 Mg, > m^o + 50 GeV 




Mq, + Mq- > rriAo + 50 GeV {i ^ j) 


B{h S7) 


(2.1 -4.5) X 10-4 (3o-) 



Here niq. = rriq^ ^ are the squark soft SUSY-breaking mass parameters in each chiral sector, which 
are common for the three generations; Mq^ are the physical masses of the squarks, M2 is the SU{2) l 
gaugino mass, and ruh stands for the mass of the corresponding Higgs boson h = h°,H°, A'^. Due 
to the structure of the couplings in the MSSM, and taking into account that the range tan (3 < 2.5 
is not favored in the MSSM, the parameter tan/3 is fixed at a high value for both channels as 
indicated. Concerning the characteristic SUSY mass parameter MsusY, it sets the scale for the 
masses of the squarks and gauginos. Notice also that our analysis incorporates a specific choice of 
sign for (> 0) partially motivated by the data on the muon anomalous magnetic moment 4. We 
adopt Bih 57) = (2.1 — 4.5) x 10^^ as the experimentally allowed range at the 3cr level [5] and 
include in our codes the MSSM computation of the branching ratio at leading order from Ref. [36]. 
In addition, we make sure that the sign of the MSSM amplitude for b — S7 and the purely-SM one 
do coincide [37]. By enforcing these experimental and theoretical constraints in our calculation, 
we automatically eschew regions of the MSSM parameter space which would artificial enhance the 
predicted qq' rates at the LHC. 

Concerning the numerical computation of the direct hs production pp bs (mainly from gluon- 
gluon fusion pp(gg) bs), some technical stumbling blocks have to be overcome. On the one hand, 
it is well-known that the presence of light quarks in the final state of such type of processes entails 
large logarithmic factors, which depend on the soft scales of the problem - precisely the masses 
of these light quarks. These large logarithms are the remnants of the truly collinear divergences 
that would arise if quarks were massless, and they turn out to be related to the non-perturbative 

''The sign // > is not essential for our numerical results. Moreover, the observable g^i — 2 also depends on the 
value of some slepton masses which do not play any role in our calculation. For the SUSY effects on g'^ — 2, see e.g. 
the excellent review [35] and references therein. 



6 



regime of QCD. Upon resummation, these terms can be factored out from the computation of 
the partonic cross section, and finally reabsorbed into the definition of the parton fragmentation 
functions. Nevertheless, we do not need to address such a detailed analysis here. Instead, we 
can include an angular cut (sin^ 9 > 0.05) to circumvent such delicate coUinear scenario. Such a 
simpler strategy should be enough to attest the fact that the direct production of 6s is, by far, a 
subleading mechanism - See Section 4 for details. 

A second (and even more subtle) difficulty is caused because the hs production threshold is 
very small. In the limit of very low external momenta (which is, by the way, the situation when 
we probe partonic \fS energies near the hs threshold), we encounter one-loop expressions of the 



where e <C mi . This kind of expression, when evaluated numerically, generates pseudo-singularities 
of the kind 0{l/e'^) in some of the intermediate steps, a fact that is reflected in the associated 
integration uncertainties. Numerical instabilities have a critical impact on the overall computation, 
which involves a large number of diagrams (some of them including various "box diagrams" , i.e. 
4-point functions). We refrain from displaying here the full list of diagrams, which is very similar 
to those shown in Figs. 1-5 of Rcf. [20], after appropriately replacing the virtual quark and squark 
contributions in accordance with the new external b and s states. Particularly subtle are the 
instabilities involved in solving the Passarino-Veltman reduction formulae necessary to obtain 
the 4-point amplitudes (a cumbersome diagonalization procedure which is extremely sensitive to 
numerical niceties) . 

A thorough analysis of such instabilities has been performed, from which we conclude that we 
can handle them by: i) changing the subroutine that undertakes the Passarino-Veltman reduction 
for the 4-point amplitudes within the LoopTools framework - we use, instead, an independent 
version presented in [38], which is also implemented in the LoopTools code but only used to cross- 
check the results in the standard setup. This way the problem is smoothed, albeit not fully solved 
since the new subroutine renders lower uncertainties in such delicate regimes; and ii) introducing 
a cut over the partonic integration domain in order to elude the neighboring region of the bs 
threshold where the numerical instabilities arise. 

4 Numerical analysis of the bs production rate 

Let us first consider the Higgs boson production/decay mechanism pp — *■ ^ — *■ 6s. The main result 
of our Monte Carlo scan is shown in Fig. 2a, which displays the maximum values of the production 
cross-section o"(pp h ^ bs), Eq. (1), for the three MSSM Higgs bosons h = h", H'', A*^ at the 
LHC as a function of m^o. On the left- vertical axis of this figure, we indicate the value of the 
cross-section (in pb) and at the same time we track the number of FCNC events (per 100 fb"-*^ 
of integrated luminosity J C dt) on the right vertical axis. One can immediately see that, in the 
large tan (3 regime: i) The maximum number of events reaches a sizeable level, which lies between 
10^ and 10^ per J Cdt = 100 fb"^ in the interesting mass region 100 GeV < m^o < 200 GeV. In 
the particular case of the h° channel, the rate of ~ 10^ events extends even farther for larger m^o 
(up to around 250 GeV); ii) Actually, for this channel, there is an approximate plateau extending 
across the mass range 300 GeV < m^o < 600 GeV, where the event rate ~ 10^ is sustained; iii) 
The other two channels H*^, can reach a similar (though smaller) maximum number of events, 
but only in the strict range 100 GeV < m^o < 200 GeV beyond which the rate decreases steadily 
and unstoppably. We can compare this behavior with the SUSY-QCD case [23] - see Fig. 2b - 
which looks similar but significantly shifted upwards. For a better understanding of these results, 
in Fig. 3a we show the corresponding FCNC decay branching ratio for the MSSM parameters 
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Figure 2: Maximum contributions to o"(pp ^ h ^ bs) in Eq. (1) as a function of m^o (at fixed 
tan/3) taking into account: a) SUSY-EW interactions with minimal flavor-mixing; b) SUSY-QCD 
interactions (see Ref. [23]). The left- vertical axis provides the cross-section (in pb) and the right- 
vertical axis tracks the number of events per 100 fb"^ of integrated luminosity. 




Figure 3: a) SUSY-EW contributions to the FCNC decay branching ratio B{h hs) as a function 
of rriAo for the MSSM parameters that maximize the various (T(pp h ^ bs) (cf. Fig. 2a); b) 
Absolute value of tana for the parameters that maximize the specific channel (t(pp — s- h*^ — >^ bs) 
where the "small a^s scenario" can be realized. 
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A" 


cr(pp h ^ bs) 


0.075 pb 


0.20 pb 


0.062 pb 


events/ lOOfb-i 


7.5 X 10^ 


2.0 X 10^ 


6.2 X 10^ 


B{h bs) 


1.26 X 10-4 


8.7 X 10-^ 


1.5 X 10-4 


T{h X) 


9.5 GeV 


1.7 GeV 


11.3 GeV 


tana 


-0.30 


0.67 


-0.01 


rriq 


880 GcV 


660 GcV 


990 GcV 




1000 GeV 


970 GeV 


1000 GeV 


Ab 


-2350 GeV 


1900 GeV 


100 GeV 


At 


-1550 GeV 


-950 GeV 


2000 GeV 


B{h S7) 


2.1 X 10-4 


2.1 X 10-4 


2.1 X 10-4 



Table 1: Maximum MSSM value of cr(pp h ^ bs) (and of the number of 6s events per 100 fb"-*^) 
at the LHC, for m^o = 200 GeV and tan/3 = 50 under the assumption that the SUSY-QCD effects 
are negligible. Shown are also the corresponding values of the relevant branching ratio B{h ^ bs) 
and of the total width of the Higgs bosons (h = h°, H°, A''), together with the values of the SUSY 
parameters. The last row includes B{h ^ S7), which is seen to lie in the allowed experimental 
range. 

that maximize the various cj(pp ^ h ^ bs). In Fig. 3b, we focus on an important feature of 
the particular h^ channel, the so-called "small a^s scenario" [30]. Wc plot there the (one- loop 
corrected) value of | tanal for the parameters that maximize the cross-section a{pp ^ bs) ^. 
The significant decrease of |tana| for m^o > 300 GeF explains the relative stabilization of this 
cross-section in that region (cf. Fig. 2a). We have also included the corresponding curve in which 
the angle a is computed at the tree-level (labelled tree-Higgs). 

In Table 1 we summarize the numerical values of a{pp h ^ bs), together with the parameters 
that maximize the production rate for tan /3 = 50 at the particular point m^o = 200 GeV. We also 
provide the value of B{h bs) and r{h — X) at the maximum of the FCNC cross-section, and 
the corresponding value of i?(b S7). 

Clearly, while the H'' and A'^ channels follow a simple monotonous behavior, the particular 
hO one behaves in a complex way. Let us further elaborate on this point. At low values of 
m^o < 300 GeV, jtanaj is large and in this region the radiative corrections can make it even 
larger (Fig. 3b). In this situation, the coupling h'^bb ~ sin a/ cos /? can be enhanced and the 
dominant production subprocess is a{pp — > h'^bb). Note that, in this scenario, r(h'' bb) is also 
enhanced whereas B{h.^ — ^ bs) is suppressed (cf. Fig. 3a). The net outcome is that the increase of 
the cross-section overcomes the suppression of the branching ratio and the final result is one order of 
magnitude larger than the tree-Higgs expectation - see Fig. 2a. At large values of rrij^o > 300 GeV, 
instead, where the renormalized value of a (i.e. the effective oteff) is significantly smaller than 
the tree-level prediction [30] (cf . Fig. 3b) , the efficiency of the production subprocess pp — ^ h'^bb 
is severely hampered. The Higgs boson production is then dominated by gluon fusion at one loop 
{gg h^) and since this mechanism is controlled by the top-quark coupling h'^tt ~ cos a/ sin/3, 
it becomes insensitive to variations in small cte//- Adding this to the fact that r(h'^ bb) is 
strongly suppressed and B{h.^ — ^ bs) is correspondingly enhanced (cf. Fig. 3a), we end up with 
a regime in which the cross-section for h° — ^ bs production is the most favored one over a fairly 
sustained range of m^o. 

^Tho values of | tariQ] shown in Fig. 3b correspond to the maximization of the SUSY-EW contributions in Fig. 2a. 
The corresponding values obtained from the maximization of the SUSY-QCD contributions (Fig. 2b) are essentially 
the same, except for the random fluctuations appearing in a Monte-Carlo computation. 
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The same effect was observed for the SUSY-QCD contributions [23] , where the production rate 
is augmented some three orders of magnitude with respect to the tree-Higgs case (cf. Fig. 2b). 
The reason for this is that the SUSY-QCD effect on r(h*^ — > bs) is proportional to cos^(a — /3) (see 
e.g. Eq. (3.5) from Ref. [22]). In the tree-Higgs case, this factor goes rapidly to zero for large my^o 
because a — > /3 — 7i"/2. However, at one loop such factor becomes cos^(Q;e// — /3) and the previous 
relation does no longer hold, so there is no such suppression. Another interesting feature to note is 
that the SUSY-EW contribution to r(h° ^ bs) is roughly one order of magnitude smaller than the 
SUSY-QCD one. This property is reflected at the level of cr(pp ^ h° — bs) (Figs. 2a,b) because 
in this regime the leading Higgs boson production mechanism is the same (gluon fusion) in both 
cases. Using our analytical approximations, we can provide a simple explanation for this numerical 
difference. If, for the sake of this estimate, we take all SUSY masses of the same order, the ratio 
between the SUSY-QCD and the SUSY-EW contributions to r(h° bs) is expected to be 

-10^^ "" ^^^ y fmwV 1 cos^(ae// - /?) 
\C(W \Vts\J \mt ) tan^/? sin^ag// 

Here wc have used the effective SUSY-QCD-induced /i'' 6 s-coupling from Eq. (3.5) of Ref. [22] 
and the SUSY-EW contribution from the WF chargino loop effects, the latter being similar to 
Eq. (2) with cos a^f f replaced by sin a^f / . The prefactor ~ 10^ comes from the numerical factors 
(2/3)^ 16^ appearing in these formulae. For tan/3 = 50 and the numerically determined values of 
S23 10~^-^ ~ 0.03 and a ~ 10~^ (corresponding to large m^o > 300 GeV), wc find that (8) is 
indeed of order 10, as confirmed by comparison of plots (a) and (b) in Fig. 2. This explains nicely 
the approximate numerical relation between the SUSY-QCD and SUSY-EW effects in the "small 
aeff scenario" from simple dynamical considerations. Finally, let us mention that, for very large 
values of m^o (namely nij^o > 800 GeV) the small cte// scenario is no longer maintained, and so 
the tree-level and 1-loop values for a tend to merge (see Fig. 3). By the same token, also the two 
•^(PP ^ h° ^ bs) curves, those for a being computed at the tree-level and at 1-loop respectively 
(cf. Fig. 2), tend to approach each other in the aforementioned limit^. It is worth recalling that 
the light Higgs mass {m-^o) reaches its upper bound (771^0 < 130 GeV) for nij^o 00. Since a also 
reaches a constant value for very large rriy^o, we can understand the reason why a{pp — > h° — > bs) 
exhibits an almost flat slope in this asymptotic regime. 

Let us now evaluate the 6s event rate attained from the direct production mechanism o"(pp 
bs). In the MSSM case, we may have both SUSY-EW and SUSY-QCD effects. The number of 
diagrams is rather large and, as mentioned above, it can be inferred from those in Ref. [20] after 
appropriate replacements of the internal and external lines. Let us first concentrate upon the purely 
SUSY-EW part. In Table 2 we present, in a nutshell, the predicted values for the direct production 
of bs pairs through gluon-gluon fusion i7(pp(gg) bs) from the electrowcak charged-current effects 
(chargino and charged Higgs boson loops). The computation of o"(pp(gg) bs) is performed within 
the parameter set that optimizes the bs production rate through pp — > h'^ — > bs (central column 
of Table 1). We notice that the non-standard contributions are tiny (of order 10~^ at most) and, 
moreover, a destructive interference arises when wc add up the chargino and charged-Higgs boson 
mediated amplitudes. The latter are suppressed by small factors of m^^/rrijj, ^wf-^-^SUSY 
compared to the SM ones. In addition, the GIM mechanism involved in the SM part is much less 
severe in the down-like quark sector as compared to its effect in the up-quark sector owing to the 
presence of the large top quark mass in the latter case. 

Unfortunately, the probability to observe "direct bs events" docs not improve in the presence of 
explicit sources of supersymmetric flavor mixing. In fact, our calculation of the SUSY-QCD effects 
on the direct bs production shows that the corresponding cross-section cannot compete with the 

similar discussion, specifically for the SUSY-QCD case, can be found in Ref. [23]. 
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partial contribution 


o-(pp(gg) bs){ph) 


Charged Higgs 


9.7 X lO-** 


Chargino 


1.1 X 10-5 


SUSY-EW 


5.7 X 10"'' 


MSSM 


1.3 X 10-^ 


SM 


1.3 X 10-^ 



Table 2: Different SUSY-EW contributions to the direct bs production for the choice of parameters 
that maximize cr(pp — h° — >^ bs) (central column of Table 1). In the absence of significant SUSY- 
QCD effects, the SM and the overall MSSM contributions are coincident. 




Figure 4: Maximum contributions to cr{pp h ^ tc) - cf. Eq. (1)- as a function of m^o (at 
fixed tan/3) taking into account a) SUSY-EW interactions with minimal flavor-mixing and b) 
SUSY-QCD interactions without electroweak effects (see Ref. [23]). 



FCNC Higgs boson decay channels. The reason is twofold: i) the mass insertions arising in the 
gluino and neutralino loops produce vertices of the guise (5, Xa) bs mb{Ai, — n tan (5) (cf. Section 
2), and so proportional to m^, whereas in the tc channel the effective vertices get a factor of mt, 
ii) the particular combination Ab — fi tan /3 becomes directly and stringently constrained by the 
experimental data on B{h 37), which impose small values of ^23 or very heavy gluino/neutralino 
masses. For example, if we take the set of MSSM parameters for which the H'^-mediated SUSY- 
QCD case is optimized (see Table 1 of Ref. [23]) we find cr(pp bs) = 1.1 x 10^"^ pb, versus 0.45 pb 
from the corresponding Higgs decay mechanism. These results can be compared with the purely 
SM contribution to direct bs production, which we find it to be <t(pp — > bs)sM = 1.3 x 10"^ pb and 
entails ~ 100 events per J Cdt = 100 fb"^. We conclude that the enhancement capabilities of the 
bs events within the MSSM are dominated by the FCNC Higgs boson decay modes rather than by 
the direct FCNC production processes. 
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h 






(t(pp ^ h ^ tc) 


8.8 X 10~^ pb 


2.0 X 10"^ pb 


events/lOOfb-i 


8.8 


2.0 


B(h^ tc) 


8.5 X 10-'^ 


2.4 X 10"'' 


T{h X) 


36GeV 


39 GeV 


tana 


0.046 


-0.11 


rriq 


300 GcV 


350 GcV 




350 GeV 


350 GeV 




-675 GeV 


-1000 GeV 


At 


20 GeV 


-75 GeV 


B{h S7) 


2.9 X IQ-^ 


2.5 X 10-^ 



Table 3: Maximum SUSY-EW induced value of cr(pp h ^ tc) (and of the number of tc events 

per 100 fb"^) at the LHC, for m^o = 300 GeV and tan/3 = 50. Shown are also the corresponding 
values of the relevant branching ratio B(h bs) and of the total width of the Higgs bosons, 
together with the values of the SUSY parameters. The last row includes i?(b — >■ S7). 

5 Numerical analysis of the tc production rate 

Let us now present the results of an equivalent maximization procedure for the tc production rate, 
f^iPP —^h—^tc). In Fig. 4a, we plot the maximum cross-sections attained at different values of 
the CP-odd Higgs boson mass, m^o. For better comparison, in Fig. 4b we show the corresponding 
maximization for the SUSY-QCD case [23]. As in the previous bs analysis, the tc production 
rate becomes again favored in the large tan /? regime: this is due to the correlations between the 
Higgs boson production and its subsequent decay (see the analytical estimates for the tc channel in 
Section 2). Further numerical details of the optimal MSSM configuration are quoted in Table. 3. 
We find that the predicted rates are rather inconspicuous (of the order of cr ~ 10~^ pb at most), and 
hence difficult to detect. We emphasize that this result is essentially triggered by the charginos. 
Neutralinos, again, have an even lesser impact. Indeed, if we recall the analytical estimates in 
Section 2 and the results from Tables 1 and 3, we observe that the optimal value of At in the latter 
(which is responsible for the neutralino contribution to the tc channel) is much smaller than the 
optimal value of A^ in the former. 

It is instructive to trace the main differences between the SUSY-EW effects on the two channels 
{tc and 6s) by using simple qualitative arguments based on the dynamical features of the MSSM. 
Fortunately, this can be immediately done from the analytical estimates presented in Section 2 
and the optimal parameter sets in each case, which can be extracted from Tables 1 and 3. We 
note that the values of tan a and tan P are quite similar in both cases. Therefore, the production 
cross-section obtained from the FCNC decays of the MSSM Higgs bosons must render essentially 
the same result and cancels in the ratio. Taking advantage of this fact, the ratio of the two cross- 
sections for h = H^, A^ (the two states available for both bs and tc final states) is essentially given 
by that of the branching ratios and also by that of the corresponding partial widths, 

(7(pp -^h^bs) B{h bs) _ T{h bs) 
(t(pp ^ h^tc) ^ B{h ^ tc) ~ r{h ^ tc) ' 

In the specific case of the heavy CP-even neutral Higgs boson, the approximate analytical form 
is obtained from Eqs.(2-5). Plugging the MSSM parameters corresponding to the optimal config- 
uration for each of the channels, we realize that the dominant contribution to bs comes from the 
chargino-mediated WF corrections, while the vertex corrections drive the main part of the tc one. 
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Therefore, we can approximate the ratio of the overall cross-sections of Eq. (9) as follows: 



which is in good agreement with the 3 orders of magnitude that arise from the exact numerical 
computation, once wc plug the corresponding values of the quark masses and the trilinear couplings. 
At from Table 1 and from Table 3. Therefore, we confirm that, in the optimal scenario, the 
difference between both channels can be well accounted for by the quark mass insertions squared 
times the ratio squared of the trilinear couplings (top quark mass and At, in the bs channel, versus 
bottom quark mass and Ai,. in the tc channel). 

In comparison, the situation for the direct SUSY flavor-changing production of tc pairs is 
remarkably distinct. As we have argued in Section 4, the enhancement factor of the relevant 
couplings is much larger in this case, being now proportional to the top quark mass mt- In 
addition, the b — ^ S7 constraints can be more easily eluded (just by an appropriate choice of the 
MSSM parameters, see Refs. [17, 20]). Numerically speaking, the direct SUSY-EW production 
of tc pairs furnishes a ~ 0.01 pb, or ~ 10^ events per J Cdt = 100 fb"^ - see Ref. [20] for a 
comprehensive discussion including the direct SUSY-QCD effects. Such optimal scenarios are 
favored by relatively low supersymmetric mass scales, namely MsuSY ~ 250 GeV, Mi,M2 ~ 
100 GeV (for the electroweak soft gaugino masses), together with large flavor-mixing values of 523- 
The compliance with the b — > S7 constraints is achieved by balancing the different SUSY-EW 
pieces involved in this process, i.e. assuming relatively light masses for charginos, neutralinos, 
charged Higgs bosons and the top squark, and also moderate values of tan (3. The dominance of 
the direct tc production mechanism holds also in the case of SUSY-QCD contributions [20]. In the 
most favorable circumstances, one can reach a ~ 1 pb from direct production, whereas the optimal 
Higgs boson-mediated output lies around a ~ 10~^ pb (cf. Fig. 4b). 

6 Discussion and conclusions 

A number of studies have been devoted to the analysis of the FCNC signatures carried by electri- 
cally neutral pairs of heavy quarks of different flavors as an strategy to unravel hints of New Physics 
in the forthcoming LHC data. These events are very rare in the SM and, therefore, their observa- 
tion could be highly revealing. Some of these studies, including ours, have focused on the possibility 
that the underlying new physics could be Supersymmetry, in particular the unconstrained MSSM. 
In this letter, we have discussed the SUSY effects on the production and subsequent FCNC decay 
of the neutral MSSM Higgs bosons {h = h°, H°, AP) into heavy quark pairs qq' = bs, tc at the LHC, 
i.e. pp ^ /i ^ qq' , and we have compared the results with the direct FCNC production mechanism 
pp — > qq'. Furthermore, we have computed the SUSY-EW corrections to (T(pp — > /i — > qq') and 
also the SUSY-QCD and SUSY-EW contributions to a{pp — > bs). These results extend the analy- 
ses previously presented for the SUSY-QCD effects on a{pp h ^ qq') [23] and the SUSY-QCD 
and SUSY-EW ones on cr(pp tc) [16-18,20]. At the end of the day, we have nicely completed 
the map of MSSM predictions for the heavy qq'-paiis produced at the LHC. 

The upshot of this lengthy investigation, which includes also the analyses from the previous 
works [17, 20, 22, 23]), is summarized in Table 4. These numerical results are obtained in full 
consistency with the stringent experimental constraints from b ^ S7. The most favorable channels 
turn out to be the following: 1) the Higgs boson FCNC decays into bs, and 2) the direct production 
of tc pairs, both of them with maximal cross-sections of ~ Ipb and dominated by SUSY-QCD 
effects. In that table, we also collect the results cr(pp 6s)max ~ 10~^pb from SUSY-QCD 
and ~ 10~^pb from SUSY-EW. Clearly, the direct production of bs pairs in the MSSM is highly 
inefficient as compared to the rate that could originate from FCNC decays of the MSSM Higgs 




(10) 
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FCNC mechanism channel 


t - c 


b-s 


Higgs decay-mediated mechanism 


SUSY-QCD 


~ 10-^ pb 


~ 1 pb 


SUSY-EW 


~ 10-"^ pb 


~ 10"^ pb 


Direct FCNC production mechanism 


SUSY-QCD dominance 


~ 1 pb 


~ 10"^ pb 


SUSY-EW dominance 


~ 10"^ pb 


~ 10"'' pb 


SM 


~ 10-7 pb 


~ 10"^ pb 



Table 4: Summary of optimal SUSY contributions for tc and bs events at the LHC, within order 
of magnitude. For the direct production, we follow the same procedure as in [20], viz. wc explore 
scenarios where there is a dominant SUSY component (SUSY-EW or SUSY-QCD) and allow a 
smaller contribution from the other. In the last row, we quote the SM prediction. 



bosons. Finally, if the SUSY-QCD part is negligible (e.g. because the gluinos are very heavy 

or the flavor- mixing terms arc too small), the SUSY-EW loops alone yield a small tc rate from 
Higgs boson decays, cT(pp h ^ te)inax ^ IQ-'^pb. Nonetheless, in this case, we also have 
<7(pp — /i — 6s)max ~ O.lpb, implying ~ 10^ 6s pairs per 100 fb"-*^ of integrated luminosity. In 
other words, in the absence of strong supersymmetric interactions, the tc signature disappears for 
all practical purposes but wc could still count on a fairly large amount of bs events triggered by 
electroweak supersymmetric sources. 

Despite the predicted number of FCNC qq' events is sizeable in some cases, it is far from 
obvious that they could be efficiently disentangled from the underlying background of QCD jets 
where they would be immersed, even in the most favorable conditions. For example, it is well 
known that the simple two-body decay /i ^ bb is virtually impossible to isolate, due to the 
huge irreducible QCD background from bb dijets. This led, long time ago, to complement the 
search with many other channels, particularly with the radiative decay /i — ^ 77, which has been 
identified as an excellent signature in the appropriate range [10]. Similarly, the FCNC Higgs 
boson decay channels may help to complement the general Higgs boson search strategics, mainly 
because the FCNC processes should be essentially free of QCD background. Notwithstanding, 
other difficulties can appear related to the misidentification of jets. For instance, for the bs final 
states misidentification of b quarks as c quarks in cs-production from charged currents may obscure 
the possibility that the 6s events can be really attributed to Higgs boson FCNC decays. This also 
applies to the tc final states [39], where misidentification of b quarks as c quarks in e.g. tb- 
production, might be a source of background to the tc events, although in this case the clear- 
cut top quark signature should be much more helpful (specially after performing a study of the 
distribution of the signal versus the background). These studies, however, go beyond the scope of 
the analysis presented in this work. 

To summarize, in this letter we have completed the calculations necessary to account for the 
FCNC-triggered qq' events that could emerge at the LHC from supersymmetric sources in the 
unconstrained MSSM. While admitting that the practical detection of these events can be difficult, 
the message from the theoretical side seems now crystal-clear: the MSSM has the potential to 
provide large amounts of heavy quark pairs from genuine supersymmetric FCNC interactions, to 
wit: in the 6s channel, the production and subsequent FCNC decay of a neutral Higgs boson entails 
potentially large enhancements (mainly in the large tan f3 regime) which could boost the predicted 
cross-sections up to a{pp — > /i — > 6s) ~ Ipb from SUSY-QCD, and ~ O.lpb from SUSY-EW, 
whilst the direct production is not enhanced at all with respect to the SM result. At variance with 
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this situation, the tc channel can be maximally enhanced from the direct mechanism, whereas the 
Higgs boson-mediated rate cj(pp ^ h ^ tc) proves to be negligible. From the SUSY-QCD side, we 
find the direct production cross-section a{pp — > te)inax ~ 1 pb, while from SUSY-EW we obtain 
o"max(pp tc) ~ 0.01 pb, the respective number of events being rather large: 10^ and 10^ tc pairs 
per J jCdt = 100 fb"^ . The tc signature should obviously be the preferred one for a more promising 
experimental tagging owing to the presence of the top quark in the final state. 
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